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ABSTRACT 

We introduce a methodology for analysis of multiwavelength data from X-ray 
selected BL Lac (XBL) objects detected in the TeV regime. By assuming that 
the radio-through-X-ray flux from XBLs is nonthermal synchrotron radiation 
emitted by isotropically-distributed electrons in the randomly oriented magnetic 
field of a relativistic blazar jet, we obtain the electron spectrum. This spectrum is 
then used to deduce the synchrotron self-Compton (SSC) spectrum as a function 
of the Doppler factor, magnetic field, and variability timescale. The variability 
timescale is used to infer the comoving blob radius from light travel-time argu- 
ments, leaving only two parameters. With this approach, we accurately simulate 
the synchrotron and SSC spectra of flaring XBLs in the Thomson through Klein- 
Nishina regimes. Photoabsorption by interactions with internal jet radiation and 
the intergalactic background light (IBL) is included. Doppler factors, magnetic 
fields, and absolute jet powers are obtained by fitting the HESS and Swift data 
of the recent giant TeV flare observed from PKS 2155-304. For the HESS and 
Swift data from 28 and 30 July 2006, respectively, Doppler factors > 60 and 
absolute jet powers > 10 46 ergs s -1 are required for a synchrotron/SSC model to 
give a good fit to the data, for a low intensity of the IBL and a ratio of 10 times 
more energy in hadrons than nonthermal electrons. Fits are also made to a TeV 
flare observed in 2001 from Mkn 421 which require Doppler factors > 30 and jet 
powers > 10 45 erg s -1 . 

Subject headings: radiation mechanisms: nonthermal — galaxies: active — BL 
Lacertae objects: general, BL Lacertae objects: individual (PKS 2155-304, Mkn 
421) 
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INTRODUCTION 



Radiation from blazars is thought to originate from a relativistic jet, closely aligned with 
our line of sight, that is powered by a supermassive black hole at the nucleus of an active 
galaxy. Blazars exhibit strong, rapidly-varying emission throughout the electromagnetic 
spectrum from radio to 7-rays. The b roadband spectral energy distributi ons (SEDs) of 
blazars consist of two components (e.g. iMukherjee et al.l 119971 ; IWeeked 120031 ): a low energy 
component that peaks at infrared or optical energies, and a high energy component that 
peaks at MeV-GeV 7-ray energies, and often extends to very high-energy (VHE) 7-rays. 
BL Lacertae (BL Lac) objects are a subclass of blazars distinguished by their weakness 
or absence of broad emission lines in a quasar-like optical spectrum. BL Lac objects are 
often further sub-divided into those that have low energy components peaking in the optical 
(known as low frequency peaked or radio selected BL Lacs [RBLs]) and those that have low 
energy components peaking in the X-rays (high frequency peaked or X-ray selected BL Lacs 
[XBLs]). 

There are two broad classes of models which are used to describe these observations: 
leptonic models, in which the emission is primarily from relativistic electrons and positrons; 
and hadronic models, in which th e emission is p rimarily from protons and atomic nuclei. 



J" 

In the leptonic models (see, e.g., iBottcherl 120071 . for a review), the low-energy component 



is interpreted as synchrotron emission from leptons, and the high-energy component as ra- 
diation from target photons Compton up-scattered by relativistic jet leptons. The target 
photons can originate fr om the synchrotron radiation (synchrotron self-Comp ton, or SSC; 
Bloom fc Marscher 1996), as well as external sources, e.g., the broad-line region (jsikora et~aL 



1994) , a dusty torus (IBlazeiowski et al.ll20001) . and/or the accr etion disk (jDermer et al. 



Dermer fc SchlickeiserT 19931). In hadronic models 



see, e.g., 



Mannheim fc Biermann 



1992; 
1992; 



Miicke fc Protheroel l200ll ; ktovan fc Dermerl l2003h . the low-energy component is still in- 



terpreted as synchrotron emission from relativistic leptons, but the high energy component 
originates from proton synchrotron or photopion production initiated by interactions between 
protons and soft photons. The high-energy synchrotron and Compton emissions create, fol- 
lowing cascading and absorption of 7-rays by the ambient photons, a second hadronicly 
induced 7-ray emission component. 

Blazar modeling has ge nerally concentrated on leptonic processes to model simultaneous 
multi- wavelength data (e.g ., Ghisellini fc Madaulll99"6 ; Li Sz Kusunose 2000 ; Bottcher fc Chiang 
20021 : iBottcher et al.ll200a Ijoshi fc Bottcher|[2007T ). In the standard time-dependent leptonic 
blazar models, the low- and high-energy components in the SEDs are simultaneously fit by 
injecting nonthermal electrons (including positrons) into the jet and allowing the electrons 
to evolve through radiative and adiabatic cooling. 
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In this paper, we use a different approach to model the SEDs of XBLs. In § 2, we fit the 
optical/X-ray vF v spectrum, and use this spectral form to deduce the electron distribution 
in the jet assuming that this emission is nonthermal lepton synchrotron radiation. We 
then use this electron distribution to calculate the high-energy SSC component. The SSC 
flux is then precisely given by the inferred electron distribution and a small set of well- 
constrained observables. The full Compton cross section, accurate for relativistic electrons 
from the Thomson through the Klein-Nishina regime, is used in the derivation. We take 
into account 77 absorption by high energy interactions with low energy jet radiation and by 
interactions with the intergalactic background light (IBL). We then apply this formulation 



to the recently observed gi ant TeV flare in the XBL PKS 2155-304 (§ 3, lAharonian et al. 
2007al : iFoschini et aDl2007h as well as to the March 2001 TeV flare obse rved from Mkn 421 
with HEGRA and RXTE (§4, lAharonian et al.ll2002l ; iFossati et aDboQsh . Allowed values of 



mean magnetic field B and Doppler factor 5d are obtained by fitting the X/7-ray spectrum 
of PKS 2155-304 in § [3] and Mkn 421 in § 4, and used to derive the apparent isotropic jet 
luminosity. The results are discussed and summarized in § 5. 



ANALYSIS 



Blazars vary on tim escales from days to months at radio frequencies (e.g.. lBottcher et al. 



2003 



Villata et al.ll20041). and on timescales as sh ort as a few hours or less at X/ 7-ray energies 



(e.g. JTagliaferri et all2003l ; lFoschini et al.ll2006l ). Photoabsorption argument s (IMaraschi et al 
19921 ) and observations of superluminal motions (jVermeulen fc Cohenlll994j ) show that the 
emission region is moving with relativistic speeds. 

We consider a one- zone spherical blob of relativistic plasma moving with Lorentz factor 
r = (1 — f3 2 )~ l l 2 . Quantities in the observer's frame are unprimed, and quantities in the 
frame comoving with the jet blob are primed, so that the comoving volume of the blob is 
V b ' = where R' b is the comoving radius of the emitting blob. The angle that the jet 

makes with the observer's line of sight is denoted by 9 = arccos /1, and the Doppler factor 
is given by Sd = [T(l — Distinct, rapid high-energy flares observed in blazars imply 

that the emitting region is confined to a small volume with a comoving variability timescale 
t' v , limited by light travel time, given by 



(1; 



For the observer, the measured variability timescale is 



tv ^ tv,min 



:i + z)Rf b 

5 D c 



(2) 
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In the remainder of this section, we compute the SSC spectrum using the electron 
spectrum derived with the ^-approximation for synchrotron radiation (§ 12.11) . We then derive 
the SSC spectrum using a model-dependent approach to obtain the electron spectrum by 
integrating over the exact synchrotron emissivity and minimizing \ 2 (§ 12. 2p . We compare the 
^-approximation and the full expression for synchrotron in § 12.31 and the Thomson and full 
Compton expressions in § 12.41 We determine the 77 absorption opacity from jet radiation (§ 
12. 5p . and present constraints based on power available in the jet (§ 12. 6p . The fitting technique 
using the exact synchrotron expression and the full Compton cross-section is described in § 
12.71 Two versions of the Heaviside function are used: H(x) = for x < and H(x) = 1 for 
x > 0; as well as if (2; x±, x 2 ) = 1 for x 1 < x < x 2 and H(x; Xi,x 2 ) = everywhere else. 



2.1. SSC Emission in the ^-approximation for synchrotron 

The ^-approximation for the synchrotron spectrum is useful in that it allows one directly 
to obtain the electron spectrum and then calculate the predicted vF v SSC spectrum, f^ sc , 
in terms of the observed vF v synchrotron spectrum, f^ yn , where e = hh>/m e c 2 is the emitted 
photon's dimensionless energy in the observer's frame. We demonstrate this approach in 
this section. The ^-approximation is also used to calculate jet power. For detailed spectral 
modeling, however, a more accurate expression is needed (§ 12.21) . 



The ^-approximation for the synchrotron flux is 



f syn ^ D c(JtUb ^ K{i ' s) (3) 
07ra7 



le.g. 



Dermer fc Schlickeiserll2002l ). Here is the luminosity distance, c is the speed of light, 



or is the Thomson cross section, N' e {pf' s ) is the comoving electron distribution, 



e{l + z) fe r 
V °d£b V es 

is a synchrotron-emitting electron's Lorentz factor, and 

B 2 

Ub = — 

87T 

is the mean comoving magnetic-field energy density of the randomly-oriented comoving field 
with mean intensity B. In eq. (jl]), z is the redshift of the object and e# = B/B cr , where 
B cr = 4.414 x 10 13 G is the critical magnetic field. Note that the magnetic field, B, is defined 
in the comoving frame, despite being unprimed. The ^-approximation is accurate for spectral 
indices approximately 2.0 < p < 3.5 (where the electron distribution is N'J^y') oc 7'~ p ). 
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The comoving electron distribution, N' e (y s ), can be found in terms of f* vn from eq. ([3]): 



<{i s )=Vln' e {i s) r , . 

The synchrotron emissivity, n syn (e) (photons cm' 3 s _1 e -1 ), is given by 



Kyni e ') - 3 CO^U B 6 



(6) 



with Hg = UB/m e c 2 . Using this and eq. §5$), one can determine the synchrotron photon 
number density, 



T3l 

I /A ^ £H • / fj\ ~ 

"si/n^ J "'syn\ c J 3 d/2/:4 2 /4 ' 

This can be converted to a radiation energy density through the relation 

3d|/f n 3rf|(l + z) 2 /^" 



(7) 



e'm e cV^(e') = 

_B D 



r 3+2 r6 / 
u L v,min u T>* i 



where we have made use of eqs. (HI) and (T5]). 



The SSC emissivity, integrated over volume, for isotropic and homogeneous photon and 
electron distributions is given by 



44sc(4) = ^a T ef ^ de' J 



, u'{e') , JV^V) 



7 



/2 



*c(g, r e 



(9) 



where for a homogeneous distribution, e' s J' S sc( e 's) = ^b e 'sfssc( e 's) an d e ' s is t ne scattered 
photon's dimensionless energy in the blob frame. In eq. ([HD the Compton scattering kernel 
for isotropic photon and electron distributions is 



F c (q, T e ) = king + (1 + 2g)(l - g) + ^ .^'. (l - q) 
L 2 (1 + T e q) 



H g; 



47 



11 ■ 



( ]Joneslll968l ; iBlumenthal fc Gouldlll970l ). where 



4/7 



The limits on q are 



r e (i - e'ji) 
1 



— — and r e = 4e 7 . 



47 /: 



< g < 1 



which imply the limits of the integration over 7': 



1 , 
— e 
2 



Imin = ^s I 1 + W 1 + TT 



e'e' 



(10) 



(12) 



(13) 
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and 

lL x = -P^V(e' - e' s ) + i 2 H{e' s - e') . (14) 



eV. 



The upper limit, eq. dHJ, takes into account Compton up- and down-scattering. In principle, 
the maximum accelerated electron energy 7 2 can be determined from particle acceleration 
theory. 

The vF v SSC spectrum is given by 
Inserting eq. Q into eq. f|T5|) and using eqs. ([8]) and ([5]) gives 



fSSC _ "D c s^ SSCK^sJ 1 "1 r\ 



^=(§) sfe/ (16) 



where 

- fee' 



1 + z ' 
1 + z 

and 

/ (1 + *K 



5d 

Using eq. (J2]) to give an estimate of the blob's radius, 



R, 



l + z 



we obtain 



cSSC _27nef ( d L V ( 1 \ f°° ^ /"*~ Mv, " 



C 3 \ 1 + z J \tv,min$T)B J Jo e ' 3 Jy . 7' 5 



de'^ di ™>J' J < ■ (17) 



The observed SSC spectrum in eq. ( fTTl) is a function of three observables, namely redshift 
z, the observed synchrotron spectrum (f^ yn ), and the variability timescale (t V}Jn i n ), and two 
unknowns, 5® and -B. We use spectral modeling to constrain these two unknowns. Note the 
quantity t Vtmin 8^yB is nearly, but not quite, a constant due to the appearance of 5-d in the 
integrand and the limits. 
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2.2. SSC Emission with Exact Synchrotron Expression 



Crusius fe Schlickeiserl ( 119861 ) have derived an expression for the synchrotron emissivit y 



from isotropic electrons in a randomly-oriented magnetic field (see also lGhisellini et al.lll988l ). 

A VSe'e 3 B 



e' J' (e 



h 



H Kirf) R i x ) 



where e is the fundamental charge, h is Planck's constant, 

A7re'm 2 c 3 



3eBh-f' 2 ' 

PTT POO 

R(x) = - d9 sin9 dt K s/a (t), 

2 JO J x/ sinO 



and K 5 / 3 (t) is the modified Bessel function of the second kind of order 5/3. The function 
R(x) can be approximated as follows: 



log(.R) =A + A x y + A 2 y 2 + A 3 y 3 + A A y A + A 5 y* 



(19) 



where y = log(x) and the values of the coefficients are given in Tabled! These approximations 
are accurate to ~ 1% in the range 10 ~ 2 < x < 10 1 ; outside this range, the asymptotic 
expressions from lCrusius fe Schlickeiserl (119861 ). 



R(x) 



1.80842 x 1/3 x < 1 



(20) 



can be used. They are accurate to better than 5% outside the range 10 2 < x < 10 1 . 
The synchrotron flux is then given by 



syn 



4nd 2 L 



4irhd 2 L 



di N'Jtf) R(x) 



(21) 



Using this and eqs. ([8]) and ([9]) one gets the SSC emissivity, 



e sJssc( e s) 



9a T d 2 L e s 



2 j2 r°° 



4<M 



de 



, f, 



syn H mo 



J D ll b JO 

and the observed flux (using eq. [T5]). 

9 (1 + z) 2 cr T e? r°° 



c'3 



c?7 



SSC 



16 Tld^^t 2 in ,/,;, 



de 



, f, 



syn p ma 



dj 



7 c(g,r) , 



(22) 



J 



12 



(23) 
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2.3. Comparison of Exact Synchrotron Expression with ^-approximation 

In Fig. [I] we compare the exact synchrotron expression for synchrotron radiation, eq. 
( 12"I]) . with the ^-approximation, eq. (j3J). In this calculation, we use 

N' e {i) = K e exp (^f) H{i - 7i) (24) 

for the electron distribution, with K e = 10 49 , 7^ = 10 2 , 7^ = 10 3 , 2 = 0.116, and various 
values of p. We use 5r> = 100, B — 10 mG, and t v ^ min = 300 s. This figure shows that the 
approximation is quite accurate near the center of the spectrum, but loses accuracy at low 
and high frequencies. At the high frequencies, in particular, the improved accuracy of the 
full expression is needed to accurately fit the VHE 7-ray data. 



2.4. SSC in the Thomson Regime 

Here we compare results using the Thomson cross section with the full Compton cross 
section derived above. Representing the synchrotron spectrum as a monochromatic radiation 
field with comoving energy density u' syn , the SSC flux can be approximated in the Thomson 
regime, analogous to eq. ([3D, by the expression 

fSSC,T no °D nrr „,/ A fl»r// /\ /nr\ 

Je a = -^—p ca T u syn ~f T N e (~f T ) , (25) 

where 



•4= »/^- (»> 

Inserting eqs. (JED and (J5j) into this, one gets 

label festhomf €a ' = — —3-^ , (27) 

where 

€ a €B<^D , . 

tsj/n — pfe ■ 

Here, if /e^ 150 and fl yn are given from observations, and the redshift (and hence the luminosity 
distance) of the object is known, then the quantity t Vjm i n BSj) is a constant in the Thomson 
regime. Due to the appearance of 5r> in the integrals of eq. (ITBl or (1231) . this dependence 
does not strictly hold when the full Compton cross section is used, but describes the general 
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behav ior. The constraint £u,mm-B$D = constant was derived previously by iTavecchio et al. 



(1998) 



In order to treat SSC emission in the Thomson regime in the formulation of § 12. II or 
I2.2[ one need only replace eq. flTUJ) by 



2 

3 v ~ V"' 4 7 r 



F T (q T ) = - (1 - Or) H ( q T ; — , 1 ) , (29) 



where qr = e' s / (47'V). A comparison between the Compton and Thomson calculations for a 
particular model with a power-law electron distribution with an exponential cut off is shown 
in Fig. [2J As can be seen, the full Klein-Nishina expression gives substantially different 
results, and would result in significantly different parameter fits. This demonstrates that the 
full Klein-Nishina expression is necessary to do accurate spectral modeling at 7-ray energies. 



2.5. 77 Photoabsorption 

2.5.1. Exact Internal Photoabsorption 

Here we calculate the 77 absorption optical depth, r 77 , due to interactions with the 
internal synchrotron radiation field. This will, naturally, be a function of the synchrotron 
spectrum, f^v n . Absorption will modify the high energy SSC spectrum by the factor 



T 77 

and the absorbed 7-rays will be reinjected to form a second injection component of high- 
energy leptons. In this paper, we neglect the additional cascade 7-rays formed by pair 
reinjection, which is a good assumption when the absorbed energy is a small fraction of the 
total 7-ray energy. 

The photoabsorption optical depth for a 7-ray photon with energy e-\ in a radiation field 
with spectral photon density n(e', //; r') is (IGould fc Schrederlll967l ; iBrown et al.lll973l ) 



dr' / djx'il — //) 



I J2/e[(l-u,') 



dd a T7 [eV 1 (l-//)]n'(eV;0 



(30) 



For a uniform isotropic radiation field in the comoving frame, n ra d(e' , fi' '; r') ~ n ra d{e')/2, so 
that 

/■oo 

(31) 



T 77 (e[)^B! b / de' a^e'.^n'^) . 
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Inserting the form of the absorption cross-section, one gets 



73/ 2 poo 

Trrti) = / & <Jfi fao) (32) 



where s = e'e^ 



+2/? + 41nw ln(l + w ) - 4L(w ) 



<K s o) = i ^| mw o - /3o mw o - i ( 33 ) 



/? 2 = 1 - l/s , wo = (1 + A))/(l - A)), and 

/"IU0 

L(w ) = dw w- 1 ln(l + w) . (34) 



Substituting the internal synchrotron radiation field, eq. (|7j), for n' rad (e f ) in eq. (132|) one gets 

= 8 94 ^ (1+ ; ) , f , (35) 

where e' and e are related by e' = (1 + z)e/5 D (see eq. [I]). 



2.5.2. 5 -Approximation for Internal Photoabsorption 



An accurate approximation for internal 77 opacity is given by the 5-function approxi- 
mation 



<T 77 (e 



for the cross section (IZdziarski fc Lightmanlll985l ). Eq. (131 p becomes 

r 77 (e x ) 



2m e c 3 R' h S 



where 



or 



251 



a T d 2 L 



(1 + zfe, ' 

f-r (i + z) Vt4 



9m r 4 / • /) 6 



Cl/i 



(36) 



(37) 



(3? 



This expression can be used to derive asymptotes and determine the importance of internal 
77 absorption. 
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Using this expression, and the fact that the blob must be transparent to 7-rays, it is 
possible to derive a lower limit on the Doppler factor. Requiring r 77 < 1 and the assumption 
that the synchrotron flux is well- represented by a power- law of index a (f^ yn oc e a ) leads to 



5 D > 



10— 1/ 



\2-2a 



TTl e C t v m in 



&Td 2 L fSyn 



(39) 



flDondi fc Ghisellinilll995h . 



2.5.3. Photoabsorption by the IBL 



Absorption by pair production will also occur due to interactions with the IBL (mod- 
ifying the observed spectrum by a factor of e~ T ~ n ). In fact, IBL absorption is found to 
dominate the internal 77 absorption for the blazars considered here where we require large 
Doppler factors to produce hard-spectra multi -TeV emission. Formulas for absorptio n due 
to the IBL ca n be fo und in lStecker et al.l (120061 . hereafter S06), IStecker fc Scully! (120061 ). and 
Stecker et al.l (120071 ). Howe ver, it has recen tly been pointed out that their approximations 
may overestimate the IBL (IDermerl l2007al hereafter D07), which is more consistent with 



the calculations of Primack et al. 



J2005I . hereafter P05); thus we use both the S06 and D07 



formulations for r 77 from the IBL. A comparison of the opacity from the IBL formulations 
of S06, D07, and P05 for PKS 2155-304 can be found in Fig. H 



2.6. Jet Power and Constraints on B 

The synchrotron emission implies a minimum total (particle and field) jet power, for 
a given Doppler factor, variability timescale, and magnetic field. The magnetic field which 
minimizes this jet power will be called B m i n . Modeling the SSC component implies the 
departure of the magnetic field and jet power from these values. 

We use the ^-approximation to deduce the nonthermal electron spectrum from the syn- 
chrotron spectrum, from which the total nonthermal electron energy and nonthermal jet 
electron power can be derived. Adding the hadron and magnetic field energy gives a jet 
power that can be compared with the Eddington luminosity. The comoving energy in the 
magnetic field is given by 

W' B = VIU B = ^— , (40) 
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while the total comoving energy in the electrons is given by 



(41) 



7i 



Substituting the electron distribution, N' e (j' s ) from eq. ([5]) into eq. (l4"T!) and recalling eq. 



7. 



12 



e(l + z)/(5d€b) , one obtains 



rripC 



Qnd 2 L 



ca T e 2 B U Bcr^fy 2 V 1 + z 



L syn i 



(42) 



where 



1 syn 



Sd<^bi' 2 /{ 1 + z ) jsyn 

&ve Bl [/(i+z) e 1 



and User = B 2 r /(8n). Thermal protons or protons co-accelerated with electrons will also 
contribute to the total particle energy The ratio of the total particle to electron energies is 
given by £ = Wp ar /W' e1 so that the energy in all of the particles (electrons and protons) is 



W 

par 



3ird 2 L !;m e c 2 



ca T ef(l + zy/ 2 U Bcr 5^ 
In our calculations, we set £ = 10. 

The total jet power in the stationary frame (i.e., the frame of the galaxy), is given by 



7/2 ^syn 



(43) 



P j = 2itB! 2 pT 2 c^ 



(44) 



flCelotti fc Fabianlll993l : ICelotti et all 120071 ) . where W tot = W' par + W' B , and the factor of 2 
is due to the assumption that the black hole pow ers a two-sided jet. This formulation gives 
parameters that minimize black hole jet power (IDermer fc Atoyanl 120041 ). though without 
taking into account whether these parameters provide a good spectral fit. The magnetic field 



that minimizes the jet power P^is obtained by solving \dPj/deB 
Doing this gives 



im e c 2 d 2 L {\ + z) 



5/2 



1 2 /7 



L syn 



for e B = B/B cr . 



(45) 



Smaller or larger magnetic field values than B cr eB,mi n are possible, but only for more powerful 
jets. The minimum Pj becomes 

2 



P 



14 



TriCe 



/2 2 



B,min 



(3T 2 U 



14 



Bar 



TIC 



6 D Tt. 



v,mm^B,mm 



l + Z 



Bar 



(46) 



13 



which makes the ratio of the jet power to the minimum jet power 



Pi 



j, mm 



, '-2 4 /— 3/2 



(47) 



where ( B = e B /e B ,min- 

It is also useful to compare the jet power with the total luminosity, which is given by 

fc 



Hot 



2nd 2 L 

r 2 



de 



(4£ 



including the beaming factor 1/2T 2 for a two-sided jet. The radiative efficiency is given by 
Ltot/ Pj- 



2.7. Fitting Technique 

We calculate the synchrotron spectrum with eq. ( l2Ti) using a model for the comoving 
electron spectrum, and calculate the SSC spectrum with eq. (123!) taking into account internal 
77 absorption, eq. (1321) . and IBL 77 absorption. The synchrotron and SSC components are 
fit separately and iteratively in the following manner (see Fig. H]): Starting values are chosen 
for the electron distribution parameters (see § 13 . 11) . 5d and B. The synchrotron spectrum 
calculated with eq. (1211 is fit to the low energy (in this case, Swift) data points by varying 
the electron distribution but keeping <5d a nd B constant, us ing a x 2 minimization technique 



i.e., the Method of Steepest Descent; e.g., iPress et al.lll992l ). Once the electron distribution 



is obtained from this fit, the high energy data are fit to the SSC component, calculated from 
eq. ( 1231) and x 2 is calculated. Then the parameters 5d and B are varied, the synchrotron 
spectrum is again fit to the low energy data to obtain the electron distribution, and another 
X 2 is calculated from the high energy data. This process is repeated until the parameters (5v 
and B) that minimize \ 2 are found. The constants of the problem are the multiwavelength 
spectral data, the source redshift z, variability timescale t Vjmin , and the lower Lorentz factor 
7J of the electron distribution. 

In order to account for systematic errors in the VHE data, which usually dominates 
the measurement errors, we performed fits by multiplying the VHE data by a factor of 
Nsysiy /v rri id) asva where v m id is the geometric mean of the VHE data, and N sys and a sys fit 
parameters that were allowed to vary within the systematic error range. N sys accounts for 
the normalization error, and a sys accounts for the spectral index error; for no systematic 
errors, N sys = 1 and a sys = 0. 
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APPLICATION TO PKS 2155-304 



The XBL PKS 2155-304, an EGRET source jHartman et alill999h and one of the bright- 



est blazars at TeV energies, has been the subject of severa 



Zhang et all 12006 



304 is z 
d l 



Osterman et al 



0.116 flFalomo et al. 



2007 



Sakamoto et al 



multiwavelength campaigns (e.g., 
2007m . The redshift of PKS 2155- 



19931 ; ISbarufatti et al.ll2006l ). giving a luminosity distance of 
540 Mpc in a cosmology where Hq = 70 km s _1 Mpc -1 , Qm = 0.3, and Q\ = 0.7. 



In July and August of 2006, the source under went several extremely brig ht flares which 



were detected by HESS (Aharonian et al.ll2007al ) and followed up by Swift (IBenbow et al. 



20061 ; iFoschini et al.l 120071 ). The measured variability timescales of the flares were as short 



as a few minutes with HESS. Swift s observing schedule did not allow it to probe such small 
timescales; however , BeppoSAX obser vations from 1996 to 1999 show X-ray variability on 
scales of ~ 1 hour (jZhang et al.l 120021 ) . Preliminary analysis of Chandra data taken simul- 
taneously with the HESS observations show them to be strongly correlated, and thus have 
X-ray variability timescales of a few minutes as well, although ground- based observations 



show the optical to be uncorrelated with the X-rays and VHE 7-rays (ICostamante et al. 



20071 ). Rapid variability does not seem to be intrinsic to only PKS 2155-304, as variability 
on timescales of a few minutes at TeV energies has been observed from Mkn 501 as well 
( Albert et al.l 120071 ) . The lack of optical correlation makes it unlikely that the optical emis- 
sion came from the blob where the flare originated, and hence we consider the Swift optical 
data an upper limit. 

The data used in our analysis are not simultaneous; the Swift ultraviolet/optical tele- 
scope (UVOT) and X-ray telescope (XRT) data are from 30 July, when it first observed 
the flares, while the HESS data are from 28 July. Detailed simultaneous HESS data have 
not yet been made available, although PKS 2155-304 was detected by HESS on 30 July 
simultaneously with the Swift detection (IFoschini et al.ll2007l ). 



Using eqs. (I2"i~j) and (12"31) . and taking taking into account 77 absorption by internal jet 
radiation (eq. [32]) and the IBL, we simulate the broadband emission during one of these 
extremely br ight flares. We fit th e fully reduced Swift and HESS data, with the Swift data 
corrected by IFoschini et al.l (120071 ) for a Galactic column density of Nh 



1.36 x 10 20 cm" 2 . 



3.1. The Electron distribution 



The electron distribution is assumed to be of the form 

-p 



7 



H h'break ~ V) + 



(49) 
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y \-(p+i) 



l' b rea k 



H{i - i break ) 



As long as j[ < i min (eq. [13]) and 72 ^> 7^, the values of -y[ and 7 2 do not affect our fits. 
This leaves us with four fitting parameters to define the electron spectrum, K e , p, and 7b reafe - 
This spectrum is physically motivated, in that radiative cooling is expected to modify the 
electron power-law index by 1. For all of out fits to PKS 2155-304, we found p = 2.7. 

The electron distribution fitting parameters are not unique, and several sets of parame- 
ters give equally good fits to a given synchrotron spectrum. The fit of the SSC spectrum to 
the 7-ray data and the parameters <5d and t v ^ min derived in the fitting routine are, however, 
essentially independent of electron distribution parameters, as long as they provide a good 
fit to the synchrotron component. 



3.2. The SSC Spectrum 



The Swift and HESS data from PKS 2155-304 were fit using the technique described in 
§ 12.71 Results of the fits can be seen in Table [2] and in Figs. [5] and El The parameters t„,mm 
and 7^ were kept constant during the fits, and the angle along the line of sight was assumed 
to be small, so that r ~ 5n for the purpose of calculating the jet power. When we did vary 
j[, we found that this had no significant effect on our fits, as long as it was sufficiently low. 
All the models presente d here have 7^ = 10 3 . Based on the systematic errors reported in 
Aharonian et al.l (j2007al ). N sys was allowed to vary between 0.8 and 1.2, and a sys was allowed 



to vary between —0.1 and +0.1. This leads to the error bars on 5^ and B in Tabled which 
are systematic error bars. 

The fit parameters are strongly dependent on t v ^ min . As t„,mm increases, the best fit Sj) 
decreases. This is mainly due to the relationship defining the SSC spectrum in eq. ( TT7|) . which 
remains roughly constant provided that the quantity t^m^S^B remains constant, which is 
true for scattering in the Thomson regime, as shown below. 



(50) 



Based on the Swift and HESS data, eq. (139]) gives the constraint 



Sd >48 



300 sec 



0.175 



Thus, the short variability in PKS 2155-304 (~ 300 sec. lAharonian et al.ll2007al ) limits the 
Doppler factor to quite high values (see also Begelman et al. 20081 ). Another constraint can 
be found from eq. (1281) by noting that the peak energy for the SSC component must be > 0.2 
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TeV (e, > 4 x 10 5 ). This leads to 



B 

Tg 



(51) 



( ITavecchio et al.lll998l : Bednarek fc Protherodll997l . Il999f ). 



The best fit magnetic field values implied by spectral modeling are a few to ~ 100 
niG, which represent a small fraction of the magnetic field, eq. (145]) . that minimizes the jet 
power. The absolute jet powers for a two-sided jet derived from these fits are ~ 10 46 - 10 47 
ergs s _1 , representing a very large fraction of the Eddington luminosity for a 1O 9 M black 
hole [LiEdd = 1-3 x 10 47 erg s" 1 ). Values this large would cast doubt on the underlying 
one-zone synchrotron/SSC and IBL model. Also note that the small values of (b indicates 
that Wp ar 3> W' B \ thus, Pj oc £, and if £ is any larger, the jet power will be larger by the 
same factor. Because the jet power is strongly particle-dominated, it is very dependent on 
the value of 7^, which is not strongly constrained by our models. However, we take = 10 3 , 
a rather large value. It is unlikely to be larger, and thus, the jet power is unlikely to be lower, 
although this part of the electron distribution is not well-constrained by observations due to 
the lack of correlated optical variability. Our choice of electron spectral index assumes the 
blob is primarily cooled by synchrotron or SSC emission, which may not be the case (see 
below) . 

Fits were performed with the IBL models of (in order of decreasing intensity) S06, 
D07, and P05. The lower-intensity IBLs do lower the Doppler factors and jet powers con- 
siderably. This may be an arguement for the lower IBLs , which has also been suggested 
from obser yations of TeV observat ions of 1ES 1101-232 (jAharonian et al.l 120061 ) and 1ES 
0229+200 (jAharonian et al.ll2007bl ). However, the jet powers are still quite large compared 
to the probable Eddington limit. The lowest Doppler factor we were able to obtain, with the 
lowest IBL and longest t v ^ min was 5-q = 58. This Doppler factor is still considerably larger 
than values measured with the VLBA from PKS 2155-304 on the parsec scale (fiobs ~ 4; 



Piner &: Edwards! 120041 ) . The 7-ray emission may, however, be formed on size scales that are 



tens to hundreds of times smaller. Note that the models differ significantly in the frequency 
region that GLAST will observe (i.e., ~ 10 MeV— 100 GeV); thus, GLAST will be quite 
useful for distinguishing between models. 



3.3. Temporal Variability 

Temporal variability can result from acceleration, adiabatic expansion, and radiative 
cooling. Nonthermal electrons that cause blazar flares are thought to be accelerated to 
high energies by a Fermi acceleration mechanism, possibly caused by internal shocks in the 



-17- 



collision of irregularities in the relativistic jetted wind. The particle acceleration timescale 
in the comoving frame for a Fermi mechanism is 



t. 



N X 

1 y a 



(52) 



where N a > 1 is the number of gyrations an electron makes while doubling its energy, and 

eB 



2iim e c 



is the Larmor frequency (e.g., iGaisserl 119901 ) . 



As the electrons radiate synchrotron and Compton-scattered radiation, they lose energy. 
The electron energy-loss rate, or cooling rate, in the comoving frame from synchrotron 
radiation is 



-7' 



syn 



-OU T Ub1 



(53) 



The SSC cooling rate using the full Klein- Nishina cross section (jJoneslll968l ; iBottcher et al. 
1997h is 



where 
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and u'(e') is the synchrotron spectral energy density, given by eq. (jSJ). With these expressions 
one can determine the cooling timescales in the observer's frame, given by 



t. 



cool 



1 + Z 

6» ' 



cool 



1 + Z 



(56) 



The acceleration and cooling timescales can be seen in Fig. [7] for Model 3 with t v , m i n = 
300 s and iV a = 10 and 1000; results for other models are similar. The SSC cooling timescale 
starts to deviate from the Thomson regime behavior at 7' > 2 x 10 5 due to Klein-Nishina 
effects, which is not seen on this plot. Overplotted on this graph is the electron spectrum 
used to fit the PKS 2155-304 data, which cuts off at 7' = 2 x 10 5 . In PKS 2155-304, the 
cutoff in the synchrotron spectrum at v m 10 18 Hz means that electrons are not accelerated 
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to extremely large Lorentz factors. If the cutoff in the synchrotron spectrum and therefore 
the electron spectrum is attributed to cooling that prevents acceleration to higher energies, 
then N a must be very large, > 10 5 . This could occur if there are more severe limitations to 
electron acceleration than implied by the simple expression given by eq. (|52|) . See § E]for 
further discussion. 



4. APPLICATION TO MKN 421 



The XBL Mkn 421 — the first blazar seen at TeV energies — has been the tar get of many 



multiwavelength campaigns (e.g., Macomb et al. 1995 ; Fossati et al. 20041, 2008h . The c am- 



paign of March 2001 was one of the longest and most complete (IFossati et al.l 120081 ). in 
which it was observed by RXTE, Whipple, HEGRA, and the Mt. Hopkins 48" telescope 
nearly continuously for seven days. During this campaign on 19 March 2001, an extremely 
bright flare was observed with all three instruments. In this section, we model this flare 
with our SSC methodology. Mkn 421 has shown variability on tim escales down to ~ 1000 s 
e.g., lAharonian et al.ll2002h . This XBL has a redshift of z = 0.03 JlJlrich et alill975l : llJlrich 



19781 ) giving it a luminosity distance of di 
mentioned in § [31 



130 Mpc using the cosmological parameters 



Fits to the March 2001 flare for different t v min and IBL formulations are shown in Table 



Hand Fig. El We fit to the fully reduced data from IFossati et al.l (120081 ). The RXTE data 



was reduced assumed a Galactic column density of 1.6 x 10 cm" , although this makes 
little difference for RXTE's bandpass. The same electron distribution as for PKS 2155-304 
was used, namely eq. (I49p . except for Mkn 421, the best fit electron power-law index was 
p = 2.2. Again, it is assumed 5v ~ T. 



As with PKS 2155-304, for a constant t % 
constraint, eq. (1391 . gives 

5d>29 



i , as 5b increases, B decreases. The opacity 



10 3 



sec 



0.2 



(57) 



The S06 and D07 IBLs are nearly identical for z = 0.04. All IBLs result in a low opacity of 
intergalactic space at such low z. Absorption by the IBL plays a small, but non-negligable 
part for low- z VHE 7-ray sources. However, the difference between the IBL formulations is 
negligible, as the solutions with different IBLs are all within the margin of error. 



The model fits of IFossati et al.l (120081 ) give 5n = 100 and B = 1 G for a blob size of 



R' b = 1 x 10 14 cm which is similar to our fits, except with a larger B. It should be noted that 
their models do not fit the TeV data particularly well. VLBA observations of Mkn 421 show 
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only weak superluminal motion, with f3 b s ~ 0.1 on parsec scales (IPiner fc Edwardsll2005l ). 
As with PKS 2155-304, the jet must slow down considerably from the 7-ray emitting region 
to the parsec-scale jet. Also as before, (b is very similar in all fits. 

The cooling timescales for Model 13 for Mkn 421 can be seen in Fig. [9j Here, the Klein- 
Nishina effects on the Compton cooling timescale are visible at 7' ~ 8 x 10 4 , indicating 
that accounting for Klein-Nishina effects is quite important for modeling this flare. As with 
PKS 2155-304, the cooling timescale is longer than the variability timescale, which poses 
problems for the SSC model. 



5. DISCUSSION 

We have presented a new approach to modeling XBL radiation in a synchrotron/SSC 
framework. We obtain the electron spectrum by fitting a model nonthermal electron distri- 
bution to the low energy radio/optical/X-ray data, assuming that it arises from nonthermal 
synchrotron processes. This electron spectrum is then used to calculate the SSC component 
and to fit the 7-ray data as a function of a small set of parameters. The number of free 
parameters, namely the Doppler factor, magnetic field, and size scale of the radiating region, 
is small enough that \ 2 fits can be performed to high quality multiwavelength data. We dis- 
cuss the results given by fitting the non- simultaneous PKS 2155-304 data (keeping in mind 
that unambiguous conclusions will require simultaneous data sets), as well as simultaneous 
Mkn 421 data. 



5.1. The Rapid Flare in PKS 2155-304 



For the giant flare in PKS 2155-304, the rapid (~ 5 min) variability implies a small 
emitting region. This creates problems for the one-zone SSC model for three main reasons: 
it requires excessively large Doppler factors, jet powers, and cooling timescales. 

For fits to the full optical/UV/X-ray spectra for PKS 2155-304, we find that very large 
Doppler factors ( <5d > 60 for PKS 2155-304) are nec essary to explain high energy emission 
with the SSC mechanism (similar to the findings of iBegelman et al.l 120081 ). However, ra- 
dio measurements of superluminal motion from PKS 2155-304 ind icate that the jet Lorentz 
factor T ~ 10, at least on parsec scales (IPiner &: Edwards! 12004 ) . The cause of this dis- 
crepancy could be that the flare radiation is produced very close to the black hole, and 
that the blobs seen in the radio have slowed down before reaching parsec scales. However, 
a deceleration episode between the inner regions near the black hole and the parsec-scale 
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region wo uld likely produce observable high energy radiation distinct from the SSC compo- 



nent (e.g., iGeorganopoulos fc Kazanasll2003l ; iLevinsonl 120071 ). If all blazars have such large 



Doppler factors, the number of blazars aligned with our line of sight becomes smaller than 
the number actually observed. 

The excessive jet powers required by our fits to the full optical/UV/X-ray and 7-ray 
spectra are another problem for the one-zone SSC model. The jet powers for PKS 2155-304 
range from ~ 10 46 - 4 x 10 47 ergs s -1 (see Table [2]); these powers can be reduced by a factor 
of 10 for a hadron-free pair jet, recalling that we assumed £ = 10. The higher values exceed 
the Eddington luminosity for a ~ 10 9 M Q black hole by a factor of a few. Even if the black 
hole exceeded 1O 9 M , it is unlikely that a black hole with a jet would be so radiatively 
efficient. La rge Doppler factors i n syn chrotron /SSC modeling were also found in analyses of 
Mkn 501 by iKrawczynski et al. £53). 



The jet powers of our fits are significantly larger than t he jet power, 6 x 10 43 e r gs s ] 



found for the SSC fit to the same PKS 2155-304 data by iGhisellini & Tavecchiol (120081 ). 
Their SSC fit is closest to our Model 5, which underfit the optical Swift data as well. The 
electron jet power in Model 5 is a factor of ~ 40 times larger than the SSC model of 
Ghisellini fc Tavecchiol (120081 ). when taking into account the fact that they calculated it for 
a one-sided jet. The main reason for this is that their synchrotron component underfits the 
optical data by a larger amount (also note that they did not use the full Compton cross- 
section, which can significantly affect derived parameter values and jet powers; see § 12. 4p . 
Even though the bolometric power at optical energies is only a factor of ~ 3 greater than the 
X-ray power, the total jet power in a synchrotron/SSC model that fits the optical as well as 
the X-rays is much greater than the power to fit the X-ray and TeV radiation because of the 
many low-energy radiatively inefficient electrons required to emit the lower energy optical 
radiation in the spectrum. It is also worth reminding the reader that, as discussed in § 13. 2\ 
the low energy part of the electron distribution — and hence the jet power — is essentially 
unconstrained by optical observations. However, it is difficult to imagine how 7^ could be 
much larger than 10 3 , so that our jet powers could be considered lower limits. Our jet powers 
are also dependent on the assumption of cooling dominated by synchrotron or SSC losses. 

Another important issue for the SSC model studied here is that one would expect 
the timescale of variability at X-ray synchrotron and TeV 7-ray energies generated by the 
same electrons would be approximately equa£]. Equating eqs. (j4]) and (126]) . assuming that 



1 Actually, the variability timescale of synchrotron X-rays should be smaller than at the corresponding 
SSC 7-rays. This is because photons with a wide range of energies are scattered by electrons with a wide 
range of energies to create the SSC emission, which tends to "smear out" the variability. The synchrotron 
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the high-energy radiation is dominated by Thomson-scattered peak synchrotron photons 
with energy e pk , implies that photons with energy e s = (1 + z)e pk e syn / (5d€b) (eq. [28]) are 
produced by the same electrons that emit synchrotron radiation with energies ~ e syn . The 
results of Model 3 for PKS 2155-304 in Table [2] imply that e s = 3.3 x 10 9 e syn , taking e pk = 
10~ 4 . Electrons Compton scattering e pk photons to make 1 TeV photons would therefore 
be radiating synchrotron photons at ~ 0.3 keV with the same short variability timescale 
as measured at TeV energies. Because the SSC cooling timescale in the Thomson regime 
scales as oc 7 /_1 , the higher energy photons observed by Swift's XRT should have an even 
shorter variability timescale than observed at 1 TeV by HESS; however, Swifts limited 
observing schedule did not allow a detailed variability study. Chandra did observe the 
flare, and preliminary anal ysis seems to indicate th e X-rays are highly correlated with the 
7-rays observed by HESS ( Costamante et al. 2007h. and strong X-ray Ay-ray correlation in 



other blazar flares ha s been observed before (e.g.. iFossati et al.ll2008l ; ISambruna et al.l 12000 



Tavecchio et al.ll200ll ). but its variability timescale is still not clear. In any case, whatever 
the source of the variability is, the timescale will be limited by the size scale of the emitting 
region. 

The one-zone SSC model for the TeV flare in PKS 2155-304 implies, from Fig. [7J 
that the radiative cooling timescales are much longer than t Vjm i n . Thus, variability cannot 
be attributed to radiative cooling, but could originate from adiabatic expansion; however, 
as mentioned above, the lack of achromatic variability in PKS 2155-304 makes this seem 
unlikely. Thus, t he va riability is not consistent with the one zone synchrotron/SSC model. 
Begelman et al. J2008h use analytic estimates to show that, with the brightness and temporal 



variability observed in PKS 2155-304, T > 50 for the radiation to avoid photoabsorption. 
They also suggest that magnetic energy density must dominate the jet power in order for the 
acceleration to be efficient, which limits T < 40; thus, the SSC model cannot explain these 
flares (however, in our simulations, the jet power is dominated by particle energy density [§ 
13.2] ). They suggest the flare is caused by Compton scattering of an external radiation field 
that must be limited to below sub-mm wavelengths. Using their inferred energy density for 
the external scattering radiation, then unless the size scale of this region is < 3 x 10 16 cm, 
this radiation should be observable, which makes it more likely to originate from an accretion 
disk. 

Are there explanations that could resolve these issues with the rapid TeV flare from 
PKS 2155-304? We discuss three possibilites: external Compton scattering, relaxing the 
assumption of homogeneity, and lower IBL energy densities. 



variability, on the other hand, is due only to the electron variability, not variability in a source photon held. 
Thus, the synchrotron's variability is not washed out as much as the SSC's. 
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The SSC mechanism involving isotropic scattering is unable to decelerate a jet blob 
from the excessive Dopp ler factors found from our fits, to the Doppler factors required 
from radio observations (IPiner &: Edwards! 120051 ) . Jet deceleration would require another 



mechanism such as external Compton scattering where jet electrons upscatter photons of an 
external radiation field. This would produce flares at GeV energies. These photons could 



originate from a slower-movi ng sheath surrounding a faster jet spine (iGhiselfini et al.l 12005 



Tavecchio fc Ghisellinil I2008T). or from a n advanced portion of the jet moving at a slower 
speed ( iGeorganopoulos fc Kazanasll2003l ). The external photons would also serve as a power 
and cooling source, which could also resolve issues with excessive jet powers and cooling 
timescales. A complete analysis would require an extension of the s ynchrotron / S S C model 
to include external Compton scattering and effects of the 77 opacity (jReimerl 120071 ) from the 
scattered radiation field to account for the full broadband SED (C. D. Dermer, J. D. Finke, 
H. Krug, & M. Bottcher 2008, in preparation). In this case, the 100 GeV - TeV radiation 
would still predominantly arise from the SSC process because of strong Klein-Nishina effects 
on external Compton components from scattered optical/UV radiation, but the additional 
77 opacity would have to be considered in the analysis. 

Another possible explanation for the extreme parameters is that the homogeneous as- 
sumption for the blob is invalid. We have assumed that the synchrotron radiation, which 
provides target photons for Compton scattering, is emitted by the same nonthermal electrons, 
whereas they could originate from a more extended region. The correlated variability be- 



tween X-rays and TeV 7-rays seen in PKS 2155-304 (IFoschini et al.l 120071 ; ICostamante et al. 



20071 ) indicate however that at least the higher energy synchrotron photons are probably 
co-spatially produced with the SSC radiation. Even in a one- zone model, the electrons could 
be strongly cooled before the synchrotron photons uniformly fill the emission region. Al- 
ternately, the entire blob might not be optically thin to 77 attenuation and the observed 
SSC emission may be from a smaller region than the synchrotron emission. Simulations with 
corrections to the one-zone approximation are necessary to investigate these possibilities. 

Recent blazar observations may indic ate an IBL energy dens i ty only slightly above the 
lower limits implied from galaxy counts (lAharonian et al.l 120061 . l2007bl ). The lowest IBL 
considered here, P05, gives significantly lower jet powers and Doppler factors for the PKS 
2155-304 flare than the other models, which may be a further arguement for a lower IBL. 
The IBLs of D07 and P05 both give jet powers below the Eddington limit for a 10 9 M & black 
hole. However, the Doppler factors are still excessive. 

To determine how certain our fit from Model 5 is, we plotted the 68%, 95%, and 99% 
statistical uncertainty contours. Also plotted is the jet power as a function of the parameters 
5d and B for Model 5 in Fig. [TOj The confidence contours trace a region that is roughly 
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constant in terms of jet power; inside the 68% contour, the power varies between a relatively 
small amount (between roughly 2 x 10 46 and 10 47 erg s -1 ). Also of note is that a relatively 
large range of Doppler factors (dp ~ 80 - 170 with 68% confidence) will give a satisfactory 
fit. Even the lowest Doppler factors allowed are quite high. The Doppler factors are too high 
for the eq. (|50|) constraint to play a part. The lines of constant jet power follow a power law 
based on B oc d^ 1 ' 5 , as one can see from eq. fpf4l) . 

p s oc R' b -%w; ar oc 5 D w; ar , 

if Sd ~ r and W' e 3> W' B . Most of the electron energy is located in the part of the electron 
spectrum below the break. Given a vF v spectral index of a\ (f e oc e ai ), and performing the 

integral in eq. (l4"3l . W par = B ai ~ 2 5^~ 4 . For constant Pj, this gives B oc S^~ ai . For our fits 
to PKS 2155-304, a± = 0.15, and we recover B oc S^ 1 ' 5 . Statistical uncertainty contours 
follow B oc S^ 2A , rather than the expected B oc 5^ (§ 12. 4p . This is due to the fact that 
scattering occurs in the Klein- Nishina regime. For 7V > 1/4, to zeroth order the scattering 
cross-section goes as approximately a kn oc e _1 . At the peak energy, eq. (??), using the 
Klein-Nishina cross section and eq. (|28[) becomes 

\fsyn( pk \1 2 \fsyn( pk \] 2 
rSSC V^synJi [J e K^syrJj 

^ B^etn B^6T ' 

which leads to B oc S^ 2 ' 6 . This is very close to the B oc ^ 2 ' 4 , considering the accuracy of 
the approximation. The internal 77 opacity, eq. (150"]) . does not significantly constrain the 
Doppler factor in Fig. (TTJ 



5.2. The March 2001 Flare in Mkn 421 



For Mkn 421, we find 5d > 30 for our model fits. VLBA o bservations of Mkn 421 show 
only weak superluminal motion, with T ~ 2 on parsec scales (IPiner fe Edwards! 120051 ). Jet 
powers are one the order of Pj > 3 x 10 45 erg s _1 , and for all of our models they do not exceed 
the Eddington luminosity for a 10 9 M black hole. The pa rameters for t his fla re in Mkn 
421 are thus much more reasonable than for PKS 2155-304. iFossati et al.l (120081 ) found the 
X-rays and VHE 7-rays to be highly correlated, indicating that they are likely emitted from 
the same part of the jet, although a highly variable external X-ray source that is Compton 
scattered by a blob would also explain this. As seen in Fig. [HI the cooling timescale is 
significantly longer than the variability timescale, as with the PKS 2155-304 flare. Thus the 
variability discussion in the previous section applies to Mkn 421 as well, and it is possible that 
Mkn 421's variability is dominated by some other mechanism. Analysis of hardness-intensity 
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diagr ams could indicate the effects of additional radiation mechanisms or variability sources 
(e.g., iKataoka et aPhoOol : ILTfe Kusunosel I2OO0I : bottcher fc ChianeJ liooj l. For example, if 
adiabatic expansion dominates the energy-loss rate of electrons with sufficiently low Lorentz 
factors, then the corresponding achromatic synchrotron spectral variability should allow one 
to distinguish this from frequency- dependent radiative cooling effects. 



5.3. Predictions for GLAST 

Fig. [TTI shows the blazar vF v flux that would be significantly detected with GLAST in 
the scanning mode as a function of observing time (see Appendix [X]). The predicted fluxes 
of PKS 2155-304 shown in Figs. E ~E have a vF v flux at 1 GeV of « 8 x 1(T 10 ergs cm" 2 
s _1 and a vF v spectral index a = 0.5. From Fig. [11], we see that GLAST will significantly 
detect PKS 2155-304 at this flux level in less than one or two ksec when integrating above 
100 MeV, and < 10 ksec when integrating above 1 GeV. Significant detection of Mkn 421 
will be achieved with GLAST in < 30 ksec whether integrating above 100 MeV or 1 GeV, 
as can be seen from Fig. [HJ For PKS 2155-304, the predicted vF v flux is quite sensitive to 
the assumed IBL (see Fig. EJ), so if the synchrotron/SSC model is valid, then the different 
IBLs could in principle be distinguished. 



5.4. Summary 

We have modeled flares in PKS 2155-304 and Mkn 421 with a synchrotron/SSC model. 
Due to the high Doppler factors, jet powers, and radiative cooling timescales, we find that it 
is unlikely that SSC emission alone can explain the giant TeV flares in the blazar PKS 2155- 
304, at least for the one-zone approximation. Lowering the IBL energy density lowers these 
quantities considerably, but not enough to change our conclusions. Although one-zone SSC 
modeling of the March 2001 Mkn 421 flare gives reasonable jet powers and Doppler factors, 
the long cooling timescales cause problems for this flare as well. Compton scattering of 
photons from an external radiation source, for example, from radiation produced in different 
regions of the jet, or loosening the one-zone approximation, might remedy these problems. 
The addition of external scattered radiation to this analysis technique will be the subject of 
future work. 
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A. GLAST Sensitivity to Blazar Flares 

GLAST Large Area Telescope (LAT) sensitivity estimates to blazar flares are updated 
using the latest GLAST LAT performance parameters^ Let u = E(GeV) represent photon 
energy in GeV. The effective area of the GLAST LAT, accurate to better than 15% for 70 
MeV < E < 200 GeV, can be written as 

a { u) = aA^ *i (A1) 

y 1 11, if 1 < u < 200, 



with Aq = 8600 cm 2 . Following the approach of iDermerl (l2007bl ). the number of source 



counts detected with the GLAST LAT above photon energy u GeV is 

V ,XAtf GeV A f 2(u a -V 2 -l)H(l-u) [max^n)]"- 1 ) 
1 U> " E GeV \ I -2a 1 - a J ' 1 ' 

Here At is the total observing time and X is the occultation factor in the scanning mode, 
7] 7 = 0.67 is an acceptance cone about the source, and E G eV = 1-6 X 10~ 3 ergs is 1 GeV in 
units of ergs, provided that the blazar uF u spectrum f e is given in units of ergs cm -2 s _1 . 
The blazar flare spectrum is assumed to be described by a single power law with vF v index 
a, so that f t = fG e vu a H(u;ui,U2), where fc e v is the vF v flux at 1 GeV, and the endpoints 
Wi < 0.07 and u 2 > 200. 

The GLAST LAT point spread function is given in terms of the angle 9(°) — 0.62u _3//4 
for 68% containment, and is accurate to better than 20% for 0.03 < u < 40. The energy- 
dependent solid angle corresponding to this point spread function is therefore AQ(u) = 
ixQ 2 = 3.7 x 10 -4 w~ 3 / 2 = ujqu^I 2 sr. The number of background counts with energy > u 
GeV is 

B(> u) = { (^Pj H(l - u) + 2 ' m ^ 1 ;"»- K+0 5 ' } , (A3) 



using the diffuse 7-ray background measured with EGRET (jSreekumar et al.l Il998l ) with 



photon index a 1 = 2.10 and coefficient k 1 = 1.37 x 10 6 ph/(cm 2 -s-sr). 



2 See LAT Instrument Performance at |http: / /www-glast.slac.stanford.edu/ 
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Fig. [IT] shows how bright the vF v flux has to be for GLAST to detect at least 5 counts 
and at n = 5 a detection, estimated through the relation n = S{> u)/ y^2B(> u) for u — 0.1 
and u — 1, corresponding to E > 100 MeV and E > 1 GeV, respectively. Results are shown 
for a flat, a = 0, and rising, a = 1/2, vF v spectrum with X = 0.2. The right-hand axis shows 
the corresponding integral photon flux for the u > 0.1, a = case in units of 10~ 8 ph(> 100 
MeV)/(cm 2 -s). The break in these curves represents a transition from a signal-dominated, 
bright flux regime where the detection sensitivity cx At to a background-dominated, dim 
flux regime where the detection sensitivity oc y/At. 
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Table 1. Coefficients for the approximation to 



Coefficient 


icr 2 < x < io° 


10° < x < 10 1 


A 


-0.35775237 


-0.35842494 


Ai 


-0.83695385 


-0.79652041 


A 2 


-1.1449608 


-1.6113032 


A 3 


-0.68137283 


0.26055213 


A 4 


-0.22754737 


-1.6979017 


A 5 


-0.031967334 


0.032955035 



Table 2. Parameters giving best fit SSC spectra to PKS 2155-304 data for different values of t var and IBL. See text 

for details. 



Model No. 


IBL 




So 


B 


Cs 


Pj 


Pj 


Ltot 


Lt 

Pi 


< 
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Imax 




reduced \ 2 






[sec] 




[mG] 




[10 46 erg s" 1 ] 




[10 42 erg s- 1 ] 




[10 15 cm] 










1 
2 
3 


S06 
S06 
S06 


30 
300 
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282 ± 22 
278±so 

168+^0 


41 ±7 
5.9±t 
2.%t\ 


0.02 
0.02 
0.03 


63 
34 
14 
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362 
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3.8 
3.9 
1.1 


1.0 x 10~ 4 

1.1 x 10~ 5 
7.6 x 10" 6 


0.23 
2.2 
14 


1.1 x 10 5 
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2.1 x 10 5 
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3000 


230+40 

124±|g 
67 ± 10 
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0.04 


3.8 
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58 ±5 
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3 x 10 40 
2 x 10 41 
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1.9 
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I 

CO 

to 

I 



Table 3. Parameters which give the best fit SSC spectra, for a given t VtJnin and IBL, to Mkn 421; parameters are 

described in the text. 



Model No. 


IBL 
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Fig. 1. — A comparison between the exact synchrotron expression (solid curves) and the 
^-approximation (dashed curves) computed for various electron spectral indices. 
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Fig. 2. — A comparison of SSC calculated with the full Compton cross section (solid curves) 
and the Thomson cross section (dashed curves) for various electron spectral indices. 
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3. — A comparison of the photoabsorption opacity for various IBL formulations for 
0.116, the redshift of PKS 2155-304. 
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Fig. 4. — Diagram demonstrating the synchrotron and SSC fitting procedure. 
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Fig. 5. — Various synchrotron and SSC fits for t v ^ m i n = 300 sec, for various IBL formulations. 
Swift UVOT and XRT data from PKS 2155-304 for 30 July 2006 (diamonds), as well as HESS 
data from 28 July 2006 (squares) are shown. The thin curves are the unabsorbed spectra, 
while the thick curves take into account 77 absorption by jet internal radiation and the IBL. 
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Fig. 6. — Similar to Fig. [5] figure compares models with different t v<m i n for the D07 IBL. 
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Fig. 7. — The cooling timescales for synchrotron (thin solid black curve) and SSC (thick 
solid black curve) as a function of electron energy for Model 5. The total cooling time from 
both of these processes is plotted as a dashed-dotted blue curve. The acceleration timescale 
is plotted as the dashed red and brown curves with N a = 10 and N a = 1000, respectively, 
and the variability timescale for this simulation (300 s) is shown as the dotted green curve. 
Also overplotted is the form of the electron spectrum, iV^Y) (solid violet curve). 
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Fig. 8. — Synchrotron SSC model fits to the March 2001 flare observed by iFossati et al. 



( 120081 ) in Mkn 421 for t VjTn i n = 10 3 sec. The triangles are limits from the Hopkins 48" 
optic al observations, the dia monds are the RXTE data, and the squares are the HEGRA 
data (jAharonian et al.ll2002l ). The thin curves are the unabsorbed spectra, while the thick 
curves take into account 77 absorption by jet internal radiation and the IBL. 
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Similar to Fig. [7] only for Model 13 of Mkn 421. 
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Fig. 10. — A plot of the jet power as a function of the parameters 5n and B for Model 5. The 
color corresponding to a specific jet power is given in the bar above, in erg s _1 . Overplotted 
are the 68% (solid curve), 95% (dotted curve) and 99% (dashed curve) confidence contours 
for Model 5. 
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Fig. 11. — The required vF v flux at 1 GeV of a high galactic latitude point source to be 
significantly detected with GLAST as a function of total time in the scanning mode. The 
thick and thin curves give the required fluxes when integrating above 100 MeV and 1 GeV, 
respectively, and the solid and dotted curves show the required fluxes for vF v spectral indices 
a = and a = 1/2, respectively. 



